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Monolayer structures made up of purely one kind of atoms are fascinating. Many kinds of honeycomb systems 
including carbon, silicon, germanium, tin, phosphorus and arsenic have been shown to be stable. However, so far 
the structures are restricted to group IV and V elements. In this letter, we systematically investigate the stability 
of monolayer structures made up of aluminium, in four different geometric configurations (planar, buckled, 
puckered and triangular), by employing density functional theory based electronic structure calculation. Our 
results on cohesive energy and phonon dispersion predict that only planar honeycomb structure made up of 
aluminium is stable. We call it ’’aluminene” according to the standard naming convention. It is a metal. Results 
of electronic band structure suggest that it may be regarded as a highly hole doped graphene. We also present 
the tight-binding model and the Dirac theory to discuss the electronic properties of aluminene. 

PACS numbers: 68.65.-k, 61.46.-w, 81.07.-b, 31.15.E-, 71.15.Mb 


Monolayer structures made up of purely one kind of atoms 
possess many novel properties different from their bulk coun¬ 
terparts. The best example is graphene, which is a hon¬ 
eycomb monolayer of carbond [2l. The other group IV 
elements such as silicon, germanium and tin also form a 
honeycomb monolayer known as silicene, germanene and 
stanene, respectively (SHU. Their geometric structures are 
buckled, where we can control the band gap by applying 
transverse electric fieldllTUTOll. They are expected to be topo¬ 
logical insulators d |5l. Very recently, silicene is experi¬ 
mentally demonstrated to act as a field-effect transistor at 
room temperature ifm . Phosphorene, honeycomb monolayer 
of phosphorus atoms, has experimentally been synthesized by 
exfoliating black phosphorus |[T2HT4ll . Its geometric structure 
is puckered. It is a direct-band gap semiconductor and pos¬ 
sesses high hole mobility ifTSll . One can also tune the band gap 
by increasing the number of layers ifTSHTTl . Another group V 
honeycomb monolayer system (buckled and puckered struc¬ 
tures) made up of arsenic atoms, called arsenene, has recently 
been predicted to be stable |[T^ . 

It is note worthy that, so far there is no report ayailable on 
monolayer structure composed of purely group III elements. 
Though there exists a study on boron based finite triangular 
structure with a hexagonal hole, it lacks translational inyariant 
symmetry to form a periodic system |[T9l . Thus, it is important 
to probe whether group III atom based monolayer systems are 
stable and then inyestigate the physical properties of the stable 
ones. 

With this motiyation, we inyestigate the stability of mono- 
layer structures made up of group III elements (B, AI, Ga and 
In) and then study the geometric, electronic and yibrational 
properties by employing density functional theory (DFT)||20l 
based calculations. We consider four different geometrical 
configurations, namely, (a) planar, (b) buckled, (c) puckered 
and (d) triangular geometries. Our DFT based phonon and co- 
hesiye energy calculations predict that among the aboye men¬ 
tioned four configurations of group III elements, only the pla¬ 
nar honeycomb monolayer made up of aluminium, as that of 
graphene, is stable. We call it ’’aluminene” in analogous to 


graphene. We also report the geometric and electronic prop¬ 
erties of aluminene with alkali metal (M = Li and Na) adsorp¬ 
tion. Further, we construct an 8-band tight-binding model and 
the Dirac theory to discuss the electronic properties of alu¬ 
minene. 

We use QUANTUM ESPRESSO packageED for perform¬ 
ing a fully self-consistent DET calculations (201 by soly- 
ing the standard Kohn-Sham equations. Eor exchange- 
correlation (XC) potential, the generalized gradient approx¬ 
imation giyen by Perdew-Burke-Emzerhof(22l has been uti¬ 
lized. We use Rappe-Rabe-Kaxiras-Joannopoulos ultrasoft 
pseudopotentials (231 for AI, Li and Na atoms which include 
the scalar-relatiyistic effect(24l. Kinetic energy cutoff of 50 
Ry has been used for electronic waye functions. We adopt 
Monkhorst-Pack scheme for k-point sampling of Brillouin 
zone integrations with 61 x 61 x I, 51 x 51 x I and 51 x 41 
X I for the triangular / buckled, planar and puckered systems, 
respectiyely. The conyergence criteria for energy in SCE cy¬ 
cles is chosen to be 10“^^ Ry. The geometric structures are 
fully optimized by minimizing the forces on indiyidual atoms 
with the criterion that the total force on each atom is below 
10“^ Ry/Bohr. In order to mimic the two-dimensional sys¬ 
tem, we employ a super cell geometry with a yacuum of about 
1 8 A in the direction perpendicular to the plane of monolayers 
so that the interaction between two adjacent unit cells in the 
periodic arrangement is negligible. To calculate the phonon 
spectra, we employ the density functional perturbation theory 
(DEPT) (25l implemented in PHonon code of the QUANTUM 
ESPRESSO package. The dynamical matrix is estimated with 
a 7 X 7 X I mesh of Q-points in the Brillouin zone. 

Stability: In order to study the stability of aluminene, we 
haye carried out phonon dispersion calculations for all the ge¬ 
ometric configurations. The fully optimized geometric struc¬ 
tures and their corresponding phonon dispersion spectra are 
giyen in Eig.I. Our calculations predict that planar aluminene 
with space group P6/mmm forms a stable structure since 
it contains only positiye frequencies for all the yibrational 
modes [Eig.I(a)]. On the other hand, the remaining geomet¬ 
ric configurations possess structural instabilities in the trans- 
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FIG. 1: Optimized geometries (top and side views) and phonon dispersion of aluminene in (a) planar, (b) buckled, (c) puckered and (d) 
triangular configurations. Planar aluminene is globally stable since the global minimum exists at the F point. There are negative phonons 
(corresponding to transverse acoustic modes) except for the planar configuration, indicating their instability. 


TABLE I: The results for optimized geometries of aluminene and 
with alkali metals (M = Li and Na) are given below. 


Structure 

Cohesive 

energy 

(eV/atom) 

Lattice Bond length (A) 
constant Al-Al Al-M 
a (A) 

M-M 

Aluminene 

-1.956 

4.486 

2.590 

- 

- 

Aluminene + Li 

-2.151 

4.517 

2.608 

2.783 

- 

Aluminene + Na 

-1.899 

4.531 

2.616 

3.197 

- 

Aluminene + 2 Li 

-2.218 

4.420 

2.552 

2.908 

2.790 

Aluminene + 2 Na 

-1.846 

4.554 

2.629 

3.217 

3.709 


verse acoustic modes due to the presence of imaginary fre¬ 
quencies (shown as negative values). Cohesive energy for pla¬ 
nar aluminene is estimated to be -1.956 eV and the negative 
value indicates that this is a bound system. Our calculations 
yield 4.486 A and 2.590 A for lattice constant and Al-Al bond 
length for the planar structure, respectively (see Table I). After 
the geometry optimization the buckled (puckered) structure 
has become AB (AA) bilayer of triangular lattice [Fig. 1(b) 
and (c)]. The relative displacement between A and B layers is 
along the (1/3, 2/3, 0) direction. Since the aluminene in trian¬ 
gular, buckled and puckered configurations are not stable, as 
observed from their phonon spectra, we do not consider them 
for further analysis. 

We have also carried out phonon dispersion calculations for 
the planar, buckled, puckered and triangular configurations 
made up of other group III elements (boron, gallium and in¬ 
dium). Our results show that there are no stable configura¬ 
tions. 

Electronic structure: We present the results of the elec¬ 
tronic band structure, density of states (DOS) and partial DOS 
(PDOS) for planar aluminene in Fig.2(a). Aluminene behaves 
as a metallic system due to the partial occupancies in the a 
as well as tt bands. We wish to note that the other honey¬ 


comb monolayer materials studied so far are either semimetal 
or semiconductor. At the high symmetric K point in Brillouin 
zone, two Dirac cones occur at energies 1.618 eV and —4.274 
eV. The former is due to the tt bonds purely made up of the Pz 
orbital and the latter is due to a bonds which has strong con¬ 
tribution from the s orbital. At the Fermi energy, there exist 
two (7 bands and we call them the cji and (72 bands as shown 
in Fig.3(a). The (72 band is purely made up of Px and Py or¬ 
bitals, while the mixing of the s orbital is not negligible in the 
(7i band. 

The electronic band structure of aluminene closely resem¬ 
bles that of graphene 1^ . but the locations of Dirac points 
in these two systems are different. In the case of graphene, 
the Dirac point lies exactly at the Fermi level where the bands 
due to the bonding orbitals (a and tt) are completely filled and 
those of anti-bonding orbitals are completely empty (a* and 
TT*). On the other hand, in aluminene, the bands due to the 
bonding orbitals (a and tt) are not completely filled and thus 
the Dirac point lies about 1.618 eV above the Fermi level. 
This is due to the fact that the number of valence electrons in 
A1 (trivalent) is smaller than that of C (tetravalent). As a result, 
aluminene can be regarded as a highly hole doped graphene. 
However, the important difference is that the Fermi energy is 
highly shifted in aluminene. Here, we wish to point out a re¬ 
markable feature of aluminene with respect to van-Hove sin¬ 
gularities. The van-Hove singularity of the a band exists very 
close to the Fermi energy at E = —0.212 eV. It may be possi¬ 
ble to tune the Fermi energy around the van-Hove singularity 
and to make the electrical conductivity very large by applying 
gate voltage. 

To understand the nature of bonding in aluminene, we have 
calculated the valence charge density [Fig.2(b)]. Except the 
hollow regions of hexagon, the charge is nearly uniformly 
extended over the aluminene plane which contributes to the 
conductivity. We have also calculated the difference between 
the charge densities of system and its atomic constituents 
[Fig.2(c)]. It clearly indicates the presence of covalent bonds 




































3 



FIG. 2: Planar aluminene. (a) Left panel: electronic band structure; 
right panel: total and partial DOS. Red curves are obtained from the 
first-principles calculation. Thin black curves are obtained from the 
tight-binding model. Circles represent van-Hove singularities, where 
the DOS is divergent, (b) Valence charge density, (c) Difference be¬ 
tween the charge density of system and their atomic charge densities 
(top and side views). 
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FIG. 3: Planar aluminene. (a) Electronic bands which cross the 
Fermi level and (b) their Fermi surfaces (curves in 2D). One Fermi 
surface is hexagonal warped, while the other two Fermi surfaces are 
of circular shape. 


between A1 atoms in aluminene. 

The transport properties of a metallic system are governed 
by the Fermi surface. The calculated Fermi surface for alu¬ 
minene is given in Fig.3(b). It consists of three closed loops, 
corresponding to the tt, cti and a 2 bands. The momentum 
vectors at which the bands cross the Fermi level along high 
symmetric /c-points are shown in Fig.3(a). Magnitudes of mo¬ 
mentum vectors along the F-K and F-M directions for the cri 
and TT bands are equal (ki = k 2 and k^ = ko). Further¬ 
more, Fermi surfaces for these two bands are circular in shape 
around the F point. Thus, the electrons in these bands behave 
as free-particles. However, the magnitudes of momentum vec¬ 
tors along the F-K and F-M directions for the (J 2 band are not 
equal {k^, 7 ^ k^). The Fermi surface for this band is highly 
hexagonal warped. We comment on these Fermi surfaces fur¬ 
ther when we analyze the Dirac theory, given below. 

Tight-binding model: Electron configuration of alu¬ 


TABLE II: The Slater-Koster parameters are obtained by fitting the 
DET based electronic band structure at the F and K points. 


Structure 

Sp 

5s 



Fssct 

Fspcr 

Aluminene 

1.618 

-1.969 

-0.815 

-1.313 

-1.596 

-1.737 

Aluminene - 1 - Li 

1.040 

-2.490 

-0.796 

-1.181 

-1.535 

-1.693 

Aluminene - 1 - Na 

0.862 

-2.256 

-0.854 

-0.985 

-1.569 

-1.620 

Aluminene -h 2 Li 

-0.128 

-3.029 

-0.750 

-1.017 

-1.601 

-1.582 

Aluminene - 1 - 2 Na 

0.089 

-3.184 

-0.867 

-1.013 

-1.395 

-1.347 

Graphene (271 

0 

-8.370 

-3.070 

-6.050 

-5.729 

-5.618 


minium is [Ne] 3s^3p^. We may construct an 8 -band tight- 
binding model including the 4 orbitals: 35 , 3py and 3pz. 
The tight-binding Hamiltonian is given by 

= ( 1 ) 

ij a 13 

where denotes the nearest-neighbor transfer integral be¬ 
tween the a orbital at the site i and the P orbital at the site j. 
We have derived the six independent Slater-Koster parameters 
(5p, 5s, Vpp^, Vppa, Vssa and Vspa) by fitting the DFT based 
band structure at the F and K points, and present them in Ta¬ 
ble IF The tight-binding model reproduces the global band 
structure well, but not precise Fermi momenta, as has been 
shown in Fig.2(a). 

Dirac theory: We proceed to construct the Dirac theory 
to describe the band structure near the Fermi energy. First, 
the Fermi surface of the (72 band is highly anisotropic [see 
Fig.3(b)]. Hexagonal warping occurs due to the hexagonal 
symmetry at the F point. It is well described by the hexagonal 
warped Dirac Hamiltonian, which was originally proposed to 
describe the Dirac fermion on the surface of the 3D topologi¬ 
cal insulator||28l, 

H = hv{kxcry - ky(Jx) + 7(^+ + (2) 

where v = 0.678a/fi is the velocity, 7 = 0.237a^ describes 
the hexagonal warping effects, Sp is the energy shift, k± = 
kx ± iky and a denotes the Pauli matrix. The energy is given 

by 

E = Sp± ^{hvkY + (n(kl - 3fc,fc2))2, (3) 

It explains well the hexagonal warped Fermi surface shown 
in Fig.3(b), where the values of the Fermi momenta k^ = 
0.269 X (27r/a) along the K direction and k^ = 0.380 x 
(27r/a) along the M direction are different (k^ > ks). 

Next, the ai band yields a circular shaped Fermi surface 
with ki = k 2 = 0.129 x (27r/a). The ai band is described by 
the Hamiltonian ^ with v = 1.996a//i and 7 = 0 . 

Third, the tt band also possesses a circular shaped Fermi 
surface with k^ = k^ = 0.375 x (27r/a). Around the F point, 
the Hamiltonian can be reduced to 

iF = (3 {dki) )Vpp7rO'x “h ^p' 


(4) 
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(b1) One Li adsorption 


(b2)Two Li adsorption 
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FIG. 4: Aluminene adsorbed with alkali metal atom (M = Li and 
Na). Geometry of (al) single and (a2) two M atoms kept at top of 
the hexagon. Electronic band structures of (bl) single and (b2) for 
two Li adsorption, and (cl) single and (c2) for two Na adsorption 
for the left panel, and the corresponding total and partial DOS for 
the right panel. Red curves are obtained from the first-principles 
calculation. Thin black curves are obtained from the tight-binding 
model, presenting good agreement for one M adsorption but not for 
two M adsorption. Enlarged version of bands and DOS for (d) two Li 
atoms and (e) two Na atoms adsorbed on aluminene near the Eermi 
level, (f) Color panel for DOS. 


We get the following energy for the above Hamiltonian, 

E = Sp ± {3 — {ak)‘^/4)VppT^, leading to a circularly shaped 
Fermi surface with the radius k = a = 

0.355 X ( 27 r/a). It is in good agreement with k = 0.375 x 
( 27 r/a) obtained from the first-principles calculations. 

Adsorption with alkali-metal atoms: In the preceding 
discussions, we have observed that aluminene behaves-like 
highly hole doped graphene and the Dirac point lies above the 
Fermi level. Then, it is expected that electron doping in alu¬ 
minene may shift the Fermi level near to the Dirac point. For 
this purpose, we adsorb alkali-metal atoms (M = Li and Na) 
on aluminene and then investigate its effects on the electronic 
properties of aluminene. Alkali-metal atoms have a tendency 
to give its valence electron due its free-electron-like behav¬ 
ior and thus they can act as an electron dopant. We consider 
three possible locations in planar aluminene, namely hollow 
(above hexagon), top (above Al atom) and bond (above Al-Al 
bond) for doping alkali metal atoms. We find that the hollow 
position gives the minimum energy configuration. Geomet¬ 


ric structures of aluminene with one and two dopant atoms in 
hollow positions are shown in Figs.4(al) and (a2). Results 
of their geometric analysis and cohesive energies are summa¬ 
rized in Table 1. 

We show the electronic band structures and DOS obtained 
from the first-principles calculation in Figs.4(bl) and 4(b2) 
for aluminene with single and two Li atoms, respectively. 
Similar results for aluminene with Na adsorption are given 
in Figs.4(cl) and (c2). We find that the Fermi level is get¬ 
ting shifted due to the adsorption of single Li and Na atoms. 
However, the amount of electron doping due to single atom 
adsorption is not sufficient to shift the Fermi level close to the 
Dirac point. We also observe an interesting feature that the 
van-Hove singularity exists at 0.116 eV for one Na adsorp¬ 
tion, which is very close to the Fermi energy. As a result of 
adsorption two alkali atoms, the Fermi level has shifted close 
to the Dirac point as shown in Figs.4(d) and (e). 

We also fix the Slater-Koster parameters for alkali-metal 
adsorbed aluminene by fitting at the F and K points (see Ta¬ 
ble II). The fitting is reasonably good for one atom doping, as 
shown in Figs.4(bl) and (cl). However, the fitting at the Fermi 
energy in the vicinity of the F point is not so good for two 
atoms doping, as shown in Figs.4(b2) and (c2). We observe 
that the bands in the vicinity of the K point move downwards 
due to alkali-metal atom doping. Consequently, Sp changes 
drastically due to the adsorption of alkali-metal atoms, while 
there are slight changes in the other parameters (see Table II). 
As a result, the adatom of alkali-metal atoms acts as the elec¬ 
tron dopant. There are only slight changes in the bands near 
the F point and they imply that the character of a bonds are 
not affected by chemical doping. 

Adding one alkali-metal atom per unit cell is not enough to 
make Sp = 0, since at the maximum it can give 1/2 electron 
to each aluminium atom. On the other hand, the Dirac cone 
exists at the K point with Sp ^ 0 when we adsorb two alkali- 
metal atoms per unit cell since they can give 1 electron to 
each aluminium atom. Unfortunately, even after the electron 
doping from the alkali metal atoms, the two sigma bands at the 
F point are not completely filled which hinder the Dirac nature 
emerging in the system. Thus, alkali metal atom adsorbed 
aluminene behaves as metal unlike the semi-metallic character 
of graphene. 

Conclusion: We have demonstrated the stability of 
graphene-like honeycomb structure made up of aluminium - 
planar aluminene - from DFT based phonon and cohesive en¬ 
ergy calculations. It is a metal with partially filled a and 
TT bands. Our band structure calculations indicate that alu¬ 
minene can be considered as a highly doped graphene. Fermi 
surface of this system contains features which correspond to 
free-electron-like and hexagonal warped surfaces. We obtain 
Slater-Koster parameters from the tight-binding model and 
construct the Dirac theory to explain the hexagonal warped 
Fermi surface of aluminene. Our investigation suggests that 
the Fermi energy can be shifted due to strong electron doping 
from the alkali-metal atoms. 

In passing we point out a remarkable feature of aluminene 
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with respect to van-Hove singularities. They emerge near the 
Fermi energy. They exist aX E = 0.856eV in the tt band and 
Sit E = —0.212 eV in the cr band. In the case of one Na atom 
adsorption the van-Hove singularity of the a band emerges at 
0.116 eV. It is intriguing to note that chiral superconductivity 
is predicted to occur at the van-Hove singularity of the tt band 
in graphene 1291 . We might expect a similar phenomenon to 
occur in aluminene. However, a study of this intriguing fea¬ 
tures is beyond the scope of the present work. 
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